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Abstract 

We propose that local quantum systems defined on expander graphs provide a simple 
microscopic model for thermalization on quantum horizons. Such systems are automat- 
ically fast scramblers and are motivated from the membrane paradigm by a conformal 
transformation to the so-called optical metric. 
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1 Introduction 



Black hole horizons are conjectured to be maximally chaotic, in the sense of maximizing 
the rate of information scrambling jUElE]. Thermalization times for sufficiently localized 
near-horizon processes are conjectured to tightly approach causality bounds, i.e. the 
information-scrambling time over a horizon patch of length L would scale logarithmically 
with the law 

r s (L)~/3 1og (f) , (1) 



where /3 is proportional to the inverse Hawking temperature and time is measured in an 
appropriately red-shifted, off-horizon frame. The arguments supporting are largely 
heuristic, based on no-cloning gedankenexperiments and ideas drawn from the 'membrane 
paradigm' of stretched horizons [U |Bj. 

In [S] we have examined these heuristic arguments to conclude that they all depend 
on a simple kinematical criterion. Namely r s in is given by the free-fall time scale 
to the stretched horizon, across the Rindler region of the near-horizon geometry. This 
means that, once is fixed, L is always bounded, so that 'fast scramblers' are necessarily 
small systems. As explained in [7J [8], the most general form of 0, applying to all kinds 
of horizons representing thermal states of hypothetical holographic dual systems, is given 
by 

T s ~ (3 log S eS , (2) 

where S e s is the specific entropy, depending on the extensive properties of the system. For 
systems smaller than a single thermal cell (the length scale set by the inverse temperature 
f3 = 1/T), the specific entropy is jus the total entropy: S e s = S. An example of such 
behavior is a near-extremal Reissner-Nordstrom black hole. On the other hand, for 'large' 
systems, S e s is defined as the entropy per thermal cell; S e s = S cc \\ = S/VT d , where VT d 
is the number of thermal cells in a volume V. A natural example would be a large black- 
brane horizon. Notice that the Schwarzschild black hole is a special borderline case, whose 
horizon size is roughly one thermal cell. 

The characterization of as 'fast' rests upon the standard intuition about informa- 
tion spread in local systems. If information is codified in some microscopic degrees of 
freedom supporting a locally conserved charge, we have transport phenomena associated 
to the conserved charge, and 'scrambling' is interpreted as the diffusive component of such 
transport. At the classical level diffusion can be described by random walk behavior of 
quasiparticle degrees of freedom, leading to mixing time scales of order 

r diS (L) = D- 1 L 2 , (3) 

where D is the diffusion constant. For systems with 0(1) dimensionless parameters, we 
can estimate that D ~ /3 is the time scale defining the coarse-graining in space and time 
over which the random walk is defined. At the quantum level, the interactions imply 
decoherence of any sufficiently localized subsystems where we may store the information. 
Hence, the classical diffusion picture should still apply provided the random-walk coarse 
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graining scale is larger than the decoherence time scale. Thus, we expect ([3]) to apply in 
order of magnitude over a wide range of conditions. 

Comparing ^ and Q at large L we see the reason for the denomination 'fast scram- 
bler'. A better characterization is obtained by focusing on the entropic form (|2]). If we 
evaluate ([3J for an extensive system with (L/(3) d thermal cells, each one contributing an 
amount of 0(1) entropy, we find 

r diff ~/3(L//3) 2 ~/3S 2 / rf , (4) 

which illustrates the 'fast' character of ^ in the limit of large entropy, as compared to 
an extensive 'slow scrambler' with the same number of degrees of freedom. 

Our discussion in this paper assumes the usual tenets of 'black hole complementarity' 
[9], namely the exterior description of the black hole can be regarded as quantum me- 
chanically complete, notwithstanding the existence of other 'complementary' observables 
which would account for the 'interior' physics. An effective model for this 'exterior de- 
scription' is assumed where Ah/^G degrees of freedom are approximately localized on the 
stretched horizon, a time-like surrogate of the event horizon with Planckian thickness. If 
the law ([T]) is applied locally on the stretched horizon, one finds that it saturates causal- 
ity [in], i.e. such a fast scrambling must operate at ballistic speed. This fact makes the 
conjectured fast-scrambling properties of horizons rather exotic and presumably require a 
highly non-local effective dynamics for the stretched horizon, despite supporting naively 
the degrees of freedom of a local system with Planck cutoff (cf. [TTJ for a recent example). 

In this note we propose a middle ground in this conceptual dilemma. Borrowing some 
standard results in graph theory we show in section 2 that there exist local scramblers 
with exponentially fast mixing, satisfying d2|). The condition being that the random- walk's 
configuration space be an expander graph (cf. p2]). We then conjecture that expander 
graphs are good local models of the stretched horizon, and we motivate this conjecture in 
section 3, using a reasoning analogous to ref. [7j. In section 4 we briefly comment on the 
general entanglement properties of quantum systems defined on expander graphs, as well 
as more intrinsic characterizations of the scrambling. Finally we spell out our conclusions 
in section 5. 



2 Local Scramblers Can Be Fast 

A precise definition of information scrambling in finite quantum systems touches upon 
many subtleties (information theoretic and dynamical alike, cf. [13J for a recent general 
discussion). For the purposes of this section, we start with a pragmatic and somewhat 
naive notion of 'local scrambler'. Namely we regard as a scrambler any quantum system 
that can be considered as a 'reservoir' for a reference 'small' system which we call the 
probe. Hence, the total Hilbert space reads 

T~L ^probc ® ^scrambler (5) 
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and we say that the probed system is a scrambler when the coupling between the two 
factors in ^ randomizes the dynamics of the probe in the limit that dim('H scram bi e r) is 
large, i.e. we have an abstract model of Caldeira-Legget type [HJ. 

We shall speak of local scramblers when a further dynamical condition is met. In 
general, the probe-scrambler coupling selects a preferred decoherence basis on "H pro be, and 
we demand that the scrambler dynamics be local in this basis. Namely, the scrambler can 
be modeled as a local many-body quantum system over some lattice, which is isomorphic 
to the decoherence basis of the probe. In this way, we model the probe dynamics as 
exercising 'quantum brownian motion' over a lattice, driven by a many-body local system 
supported on the same lattice, so that dim('Hp ro be) ~ log dim(?4crambier)- F] 



2.1 Walks 

Upon integrating out the scrambler degrees of freedom, the effective probe dynamics 
is described by a non-unitary master equation for the probe density matrix p pro be- After 
an appropriate coarse-graining in space and time, the master equation often reduces to a 
discrete Markov process over the decoherence basis of ^ pro be (cf. for example [15]). 

Let us represent a decohered probe state as a diagonal density matrix p vw = p v S vw 
in the decoherence basis, where the components p v satisfy positivity and normalization, 
Y^v=iP" = 1> t nus defining a probability distribution. Pure basis states |t>o)(^o| are given 
by sharp distributions p v = 5 VVo , whereas maximal mixing is represented by the uniform 
distribution u, defined by the components 

u v = ±. 

The Markov process is determined by a N x iV matrix M vw of transition probabilities 
between pairs of such basis states, so that on the (n + l)-th step 

p (n+l) = Mp (n) _ (g) 

Microscopic reversibility requires M to be symmetric. We can codify the matrix of tran- 
sition probabilities between basis states by means of an abstract graph consisting of iV 
vertices, one for each basis state \v), and edges (vw) between pairs of distinct vertices for 
each no n- vanishing off-diagonal matrix element M vw ^0, w^ro. 

With no essential loss of generality, we may impose the discrete analog of homogeneity, 
so that every vertex has the same degree of connectivity with the rest of the graph, a fixed 
coordination number k < N. Further imposing the analog of isotropy, we assume that 
all transition probabilities out of a given vertex are equal to one another. Finally, the 
probability of no jump from a given site is denoted r e [0, 1], and controls the diagonal 
elements M vv . All in all we have 

M = A + r , (7) 

K 

1 It is important to keep in mind that the 'lattice' of decohered states is not necessarily tied to 'real' 
space, and could arise as an effective description in, say the index space of a matrix model. 
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where A is the adjacency matrix of a simple, A;-regular graph of N vertices Gn,u- The 
adjacency matrix of such a graph is known to be symmetric with spectrum in the interval 
[—k, k\. The uniform distribution is always an eigenvector of both A and M 

Au = ku , Mu = u , (8) 

which is non-degenerate if and only if the graph is connected. The lower bound for 
eigenvalues of the adjacency matrix, —k, is attained if and only if the graph is bipartite. 
It is useful to define the Laplacian matrix of the graph 

A^l-ArM, (9) 

a positive-definite matrix with eigenvalues on the interval [0,2], which gives a discrete 
approximation to the continuum Laplacian operator —V 2 . The uniform distribution is a 
zero mode of the Laplacian, Am = 0, and the second largest eigenvalue of A defines the 
gap A, i.e. the lowest positive eigenvalue of A. 

In terms of the Laplacian matrix, the Markov process can be rewritten as a discrete 
diffusion equation 

p (n+l) _ p (n) = _(1 _ r ) A p (n) , (10) 

whose explicit solution is given by the discrete version of the heat kernel, 

pW = M> (0) = (1 - (1 - r)A) n j9 (0) , (11) 

in terms of the initial state distribution p(°\ 

Hence, all scrambling properties, at least as far as they are parametrized by the deco- 
herent (effectively classical) walk of the probe, are determined by the adjacency matrix 
A or, equivalently by the graph Laplacian matrix A. 

For our purposes, the main property of interest is the efficiency of the approach to 
the uniform distribution as a function of the size of the graph. Let us look at the norm- 
difference between the n-th step distribution and the uniform one 

\p {n) -u\ = |M> (0) - u\ = \M n {p (0) -u)\. 

Since ^2 v p^ = 1, it turns out that p^ — u is always orthogonal to u, and we can bound 

\p in) -u\ < \m\ n \p (0) -u\ < \m\ n , 

where \m\ is the largest eigenvalue (in absolute value) of the matrix M — 1. At this point, 
it will be technically convenient to work with sizable no-jump probabilities r > 1/2 which 
ensure positive-definiteness of M, so that \m\ = rh is the largest eigenvalue of M — 1, 
without regard to sign. In particular, we make the conventional choice r = 1/2, known 
in the graph theory literature as the 'lazy random walk'. We shall see further down the 
discussion that these choices are not too essential from the physical point of view. 
For the lazy random walk, m = 1 — | A and the bound becomes 

\ p (n) _ u | < (i_ i A)" . (12) 
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In other words, the rate of approach to the uniform distribution is controlled by the 
Laplacian gap A. 

Our main interest is the behavior of A as the graph becomes very large. As a simple 
example, consider the case of a regular ci-dimensional hyper-cubic lattice (Z;) d of toroidal 
topology, with degree k = 2d and vertex size N = l d . The Laplacian spectrum is 

A <f = \H - ""(WO) • 

i=l 

with S e (Zi) d and gap A = Atx 2 /I 2 = 4ir 2 /N 2 ^ d . Hence, the p^ distribution becomes 
exponentially close to uniform for n 3> n mix = N 2 ^ d . f\ Since N ~ log dim?4 cram bier 
controls the entropy of the scrambler, we recover the diffusion time scale 

T d[S ~Pn mix = PN 2 / d ~l3S 2 / d , 

as expected. 



2.2 Expanders 

Much more interesting is the case of random walks on graphs having the property 
that the gap remains bounded above zero, uniformly in the limit iV — > oo. Such graphs 
are known as expanders and feature prominently in graph theory, on account of their 
numerous applications. 

A central result of expander graph theory is precisely the exponentially fast mixing of 



random walks. In particular for our lazy random walk we can bound (12) by 1/N a with 
a > 0, by setting n > C log(iV) where C is positive numerical constant proportional to 
— a I log(l — Aqo/2) and Aoo > is the N — >■ oo limit of the gap. 

Hence, we reach the main lesson of this section: perfectly local random walks can 
behave as fast scramblers provided the graph is an expander. 

Not-so-lazy random walks with r < 1/2, hopping on A-matrices with sufficiently 
negative eigenvalues, can have transition matrices M with large negative eigenvalues which 
control the mixing properties. The extreme case of this kind is a bipartite graph, for which 
the long-time walk just oscillates between the two components, while well-randomized in 
each one of them. Therefore, the lack of precise mixing is related in this case to exact 
or approximate discrete symmetries of the graph. We can handle these cases by defining 
scrambling in terms of a single randomized component, which is equivalent to increasing 
the time coarse-graining to even time steps. 

We have defined expander graphs in terms of the spectral properties of the Laplacian. 
There is a geometrical interpretation of the Laplacian gap, as a consequence of the so- 
called Cheeger inequality 

h 2 

— < A < 2h , (13) 

2 If we demand that — u\ < 1/N a for some 0(1) positive number a, we actually need a slight 
overkill with respect to n mix , i.e. we require n > CN 2 / d log(iV) with C another 0(1) constant. 
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Figure 1: The infinite regular tree, or Bethe lattice, is the ultimate expander. Its adjacency 
matrix has a continuous spectrum in the interval [— 2\/k — 1, 2\/k — 1]. The figure shows the 
case k = 3, a tessellation of the two-dimensional hyperboloid. 



where the Cheeger constant, h, is the (normalized) edge expansion of the graph, defined 
as the formal ratio of 'area/volume', 

. # edges [OS] 
h = mm , (14) 

where the minimum is taken over all graph bi-partitions G = S U (G — S), with the 
proper subgraph S having vertex size Ns smaller than half of G, i.e. N$ < N/2. In the 
numerator we count the total number of edges connecting the two sides of the partition, 
thus measuring the boundary volume in edge units. Therefore, an expander graph family 
has a bounded finite h in the large N limit, implying that large subsets have a boundary 
volume scaling as their bulk volume. We recognize expanders as the discrete analog of 
hyperbolic geometry. 

For example, the simplest and in many ways the most perfect expander is a regular 
tessellation [TU] of the Hyperbolic space H fc_1 known as the Cayley tree or Bethe lattice 
(cf. figure 1). The metrical properties of the fc-regular Bethe lattice give an approximation 
of the (k — l)-dimensional unit hyperboloid 

ds^ h -! = dr 2 + (sinh r) 2 dVt\_ 2 , (15) 

with lattice spacing of order unityj^] We can introduce hyperbolic spaces of arbitrary 
curvature radius R, just scaling this metric by a factor of R 2 , and the corresponding 
Bethe lattice would have spacing of O(R). 



The discrete diffusion equation (10) is also reminiscent of the continuum diffusion 
equation (cf. [IS]): 

d t p = DV 2 p, (16) 



3 Continuum models of finite expander graphs can be constructed as compact hyperbolic manifolds of 
the form H fc_1 /r, with T a free discrete subgroup of isometries. 
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where D is the diffusion coefficient, p(t,x) is the continuum diffusing density and —V 2 is 
the Laplacian. For the hyperbolic space with curvature radius R it reads 



where — Vg fc _ 2 is the Laplacian on the unit sphere, R is the curvature radius of the 
hyperbolic space and we have defined z = rR. 

To investigate the speed of diffusion we look at the heat kernel 

K(z, O, t) = (z, n\ exp (t D- 1 V^) |0) 

describing the diffusion of a distribution which begins sharply localized at the origin 
z = 0. Since that distribution is spherical, we may rewrite the diffusion equation as a 
one-dimensional equation for the radial density 

f(z,t)=K(z,t)(sinh(z/R)) k ~ 2 , 

which reads 

d t f + d z (vf) = Dd 2 z f. 
This equation has the Fokker-Planck form with a drift velocity 

(fc-2) D 
V{Z) tanh(z/R) R ' 

which approaches exponentially fast the limiting value v^ = (k — 2)D/R as z > i?. 

In comparing to the discrete diffusion over an expander graph, we recall that the lattice 
spacing is or the order of the curvature radius, and assuming the free-flight velocity is of 
(9(1) both are of the order of the time step /3. Hence the continuum problem approximates 
the discrete one for D ~ R = G(f3). The drifted radial distribution is peaked around 

•^peak(^) 

~ Voo t + 0(y/tj3), with Voo of order one, implying a radial spread of ballistic type 
for the hyperbolic diffusion. 

The same ballistic character can be argued directly in the discrete random walk on an 
expander graph. Defining discrete spheres as the set of vertices at edge minimal distance 
n from a fixed central point, and balls B n as the interior of those spheres, we can reason 



as follows. Equation (14) implies that the number of edges coming out of a ball B n is 
larger than hkN n where iV„ is the vertex size of B n . Since at most k edges can end on a 
vertex, the number of vertices in the B n+ \ ball is at least N n + hN n = (1 + h)N n . Thus 
N n > (1 + h) n and we see that the vertex volume of concentric balls grows exponentially. 
This means that the size of a iV-vertex ball is of order log N, the same as the random-walk 
time for mixing on this sphere, showing that the spread of the walk is indeed ballistic. 
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3 The expander at the edge of the horizon 



Having discussed how abstractly denned local scramblers can be fast, we now come 
back to original black hole problem. An intriguing connection arises between expanders 
and black-hole horizons through the universal properties of the near-horizon geometry and 
in particular the natural description of stretched horizons in the so-called optical frame, 
obtained by removing the |<?oo| factor in any static metric by a conformal transformation 
(cf. for example [TH] for a recent review). 

3.1 Expanders and de Sitter 

The basic observation starts from the geometrical properties of expander graphs de- 
scribed in the previous section, namely the effective description of brownian motion as 
hyperbolic diffusion. An immediate connection to a gravitational system is found by not- 
ing the global conformal map between the space R x H n and the static patch of de Sitter 
space-time, dS n+ i, with metric 

ds 2 dS = (1 - u 2 ) alt 2 + (1 - u 2 Y l du 2 + u 2 alttl^ . 

The corresponding optical frame is given by 

ds 2 ~ = -alt 2 + — — + U dn 2 n _ x = -dt 2 + dr 2 + (sinhr) 2 dn 2 n _ t , 

which is globally identical to R x H n under the coordinate change sinh r = u(l — u 2 )^ 1 ^ 2 . 
Under the conformal map, null trajectories on dS n+ i are mapped to null trajectories on 
R x H n , which project to geodesic arcs on H n . Therefore, if the hyperbolic diffusion is 
realized as a random walk with piecewise geodesic arcs of unit length (in units of the 
curvature radius), it maps to a random walk on the static patch of de Sitter space-time 
made out of piecewise-null free paths. In other words, if we model scrambling on the 
static dS patch as the piecewise-null scattering of a localized probe with free path time 
scale given by the curvature, such a system behaves as a fast scrambler. 

If the hyperbolic random walk is regularized to N vertices, the whole system is scram- 
bled in a time of order log N. Since the vertex density is of 0(1) in units of the hyperbolic 
curvature, N also scales as the hyperbolic volume, which in turn scales like the boundary 
'area'. Thus, restoring units we have a time scale of order H~ x log Sgh with H the Hubble 
parameter and S'gh the Gibbons-Hawking entropy of de Sitter space-time. 

The understanding of hyperbolic geometry as the key to fast scrambling is also manifest 
in a different mechanical model for fast scrambling of dS space-time. Let us suppose that 
the probe, instead of executing a brownian motion with time scale of 0(1) in units of 
the curvature, glides freely through dS space-time except for localized scattering at the 
Planckian stretched horizon, which is regarded as a set of Sgh scatterers with Planckian 
cross section. In order to derive a bound on scrambling speed, we take the probe glides 
between horizon scatterings as null paths on dS, so that the complete trajectory maps to 
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a piece-wise geodesic trajectory on the optical frame. This defines a Hadamard billiard 
obtained by carving S GH geodesic spheres out of H n , with cross section of 0(1) in units 
of the curvature radius [21 J . The scrambling time for a probe with collision resolution of 
the order of one scatterer in such a billiard is given by the crossing time of the billiard, 
i.e. its diameter t s ~ (3 log Sgh, where we have restored the radius of curvature j3. 

We thus make contact with the chaotic scrambling models of [7], exhibiting the 
Hadamard billiard probe-scattering and the expander graph probe- diffusion as two ex- 
treme models of fast scrambling in the de Sitter static patch, depending on whether we 
confine probe interactions strictly to the Planckian stretched horizon, or we conceive of 
an 'extended' stretched horizon covering the whole dS patch, as far as probe interactions 
concerns. 



3.2 Local analysis 

The previous discussion exhibiting the global map between the de Sitter static patch 
and the hyperbolic geometry can be generalized to a local approximate map for any regular 
horizon. The basic observation is again the conformal optical map between R x H n and 
Rindler space-time, with metric 

ds^ = -y 2 dt 2 + dy 2 + dt 2 , 

where the t coordinates parametrize R n_1 . The optical frame then yields R x H n again, 
presented now in the upper-half plane coordinates 

,2 ,2 dy 2 + di 2 
ds 2 ~ = -dt 2 + . 

fun y2 

Therefore, cutoff versions of Hadamard billiards or expander graphs on the hyperbolic 
upper-half plane representation yield fast-scrambling models for Rindler space-time, which 
in turns gives a local approximate description of any regular horizon. 

To see this, notice that near a regular, static horizon we may choose ((i+2)-dimensional 
coordinates exhibiting a two-dimensional Rindler factor in the near-horizon geometry 
X cZ+2 , of the form 

dsl = g, v dx»dx» « F(p, y) dt 2 + dp 2 j + H tj (p, y) dy 1 dy> , (17) 

where the horizon is conventionally located at p = 0. The functions F(p,y) and Hij(p,y) 
are assumed to have smooth expansions around p = 0. The normalization of the time- 
time component takes into account the fact that the Hawking temperature is assumed to 
be 1//3, so that the imaginary-time periodicity t = t + ij5 is equivalent to the absence of 
conical singularity in the Euclidean section. 

The optical manifold X is defined in the exterior region p > through the conformal 
transformation 

ds\ = Q 2 (p,y)ds\, Q(p,y) = 2 ^y/F(p,y) , (18) 
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with near-horizon asymptotic form 



ds\ « -dt 2 + dz 2 + e 47rz//3 7ii (y) dj/Mj/ 3 ' , (19) 

where 7^ = (Hij/F) p=0 + O(p 2 ), and we have defined a Regge- Wheeler or 'tortoise' radial 
coordinate p = e - 27rz/ P, in terms of which the Rindler region extends roughly along the 
interval z G [0, 00], with the horizon sitting at z — 00. 

Up to terms exponentially suppressed at large z, the near-horizon optical metric is 
'asymptotically locally hyperbolic', with obvious '^-expanding' properties in the sense 
of the last section. The spatial metric d£ 2 = 7^ dy % dyi is induced at the edge of the 
Rindler region, and it is related to the physical induced metric at the horizon by a regular 
conformal transformation. Since 7^ defines a conformal compactification of the horizon 
boundary, the obvious analogy with the AdS/CFT correspondence was exploited in Ref. 
[T7] to propose a dual Euclidean CFT description on the 7-geometry. In this paper, we 



take a different route and emphasize the relation between ( 19 ) and microscopic models of 
the stretched horizon. 



3.3 The optical stretch 

In this section we explain why the optical metric is a useful tool in any discussion of 
near-horizon dynamics. In general, in the context of 'black hole complementarity' [9], we 
shall define the stretched horizon as the limit of applicability of Low- Energy-Field-Theory 
(LEFT). Any such description involves a finite number of local fields with Hartle-Hawking 
boundary conditions at p = 0, which are equivalent to the specification of a thermal state 
at locally measured temperature (/3y/— goo) -1 - Let us parametrize the effective action of 
any such LEFT in the schematic form (here for a scalar field degree of freedom \]/) 



5, 



LEFT 



X d+2 



72 , ^ or v 1 , i\/r2\ ,t, , ^ 



2 * r x + + M ) * + h> A^ om 



■ (20) 



In this expression, the Laplacian has been corrected by a conformal coupling to the 
Ricci-scalar i?[X], so that the remaining terms give the strength of conformal symmetry 
breaking. M is a typical mass, standing as a representative of whatever relevant operators 
we may have turned on, and the remaining sum accounts for characteristic irrelevant 
operators with conformal weight A > d + 2, dimensionless coupling A and scale A, playing 
the role of Wilsonian UV cutoff for the LEFTj^] Consistency of the Wilsonian procedure 
requires M < A. Since the Hartle-Hawking quantum state to be stored on the fields \l/ 
has energy scale 1/(3, the Wilsonian consistency also requires the constraint 1/(3 < A. 



At this point we are free to perform field redefinitions in (20) if they serve the purpose 
of illustrating the physics, just as it is sometimes natural to switch back and forth between 
so-called Einstein and string-frames in the graviton-dilaton system arising in low-energy 



4 We adopt the 'naturality' convention in which the strength of the coupling is controlled by A, i.e. 
A = 0(1). 
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string-theory Lagrangians. In the same vein, we shall rewrite (20) in the 'optical' frame 



defined by the conformal transformation (18), together with the field redefinition 



Using the identity 



x 4(d + l) V V 4(d + l) 



we finally obtain 
Sleft = " 



"(21) 

The result is a conformally-invariant effective action, perturbed by relevant and irrelevant 
operators with effective position-dependent dimension-full parameters M e g(z) = Q(z)M 
and A e ff(z) = fl(z)A. Since Q(z — > oo) — )■ 0, all these dimension-full parameters are scaled 
to zero in the horizon limit. Hence, the z — > oo limit is a zoom into the UV structure of the 
LEFT, isolating the conformally-invariant dynamics of whatever field degrees of freedom 
exist below the Wilsonian cutoff A, while at the same time blowing up the irrelevant 
operators heralding the approach to the UV cutoff]^] Since the thermal state of interest 
has scale 1//3, the running of the effective dimensionless coupling of irrelevant operators, 
defined at scale 1/(3, is given by 

Acff( " } = piRjra • ( 22 ) 

Since A = 0(1), the effective position-dependent coupling enters strong-coupling at the 
point za determined by 

tt(z A )A = 1/(3 . (23) 

Therefore, z = z\ marks the end of the LEFT domain of applicability, and thus defines 
the edge of the A-stretched horizon. The minimal choice for A is on the order of the 
Planck mass, Mp, defining the standard notion of a Planckian stretched horizon at the 



corresponding solution zp of (23), but other choices are possible. Consider for example 
a perturbative string theory with string mass scale M a , which defines a stringy stretched 
horizon at optical depth z s . In general we have A = M s = gfMp with g s < 1 the string 
coupling and a > 0. Solving the equations 1/(3 = Q(z s )M s = fl(zp)Mp yields the optical 
depth difference 

(3 , fMp\ (3 



5 The conformal mass term is negative, approaching — (ird//3) 2 as z — > oo. Since the X-Laplacian is 
gapped in the same amount, only the zero mode of \f r is tachyonic. This mode is seriously non-normalizablc 
in the original X-frame field ^, so it does not satisfy the Hartle-Hawking boundary conditions and should 
be treated as non-dynamical in the LEFT. 
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In this case, despite the fact that the string theory remains weakly coupled at z ~ z s , the 
effective mass of the whole Hagedorn tower of massive string states comes down below 
the 1//3 scale as z < z s . Therefore, the dynamics is expected to become non-perturbative 
because of the exponential growth of field species. 

Another interesting case is that of a 'thick' stretched horizon, arising for A ~ i.e. 
a situation where there is essentially no LEFT description at all. This corresponds to the 
bulk dynamics of holographic gauge theory duals with 't Hooft coupling of 0(1). Such 
models are akin to large-iV, weakly coupled matrix models with possibly complicated 
interactions. 



The physics of (20) in the radial domain <C z <C z\, deep within the Rindler region 
but still clear from the stretched horizon, consists of a conformal thermal state. Since 
the optical metric is ultra-static, the thermodynamic functions of such a thermal state are 
extensive on the hyperbolic spatial geometry. For example, the entropy of \l/-field states 
contained in Rindler, up to optical depth z, scales with the optical volume 

~ V z $~ d ~ X , (24) 

where 

r z p 2irzd/fi 

V Z = V dz' e Wd//3 « V , (25) 

Jo 27Td 

and V = V 1 is the volume in the jij metric. Hence, we have a uniform excitation density of 
(9(1) degrees of freedom per thermal volume f3 d+1 of optical bulk, although the 'expanding' 
nature of the geometry forces all extensive quantities to become dominated by the large 
z-cutoff (states pile up in the last layer of 0(/3) thickness). The old result of 't Hooft's 
brick-wall model [20] follows in this language as the (approximate) matching of the entropy 
in sub-Planck Hawking radiation with the Bekenstein-Hawking entropy of the black hole: 

S B R-V Zp r d ~ 1 - (26) 

Again, this entropy is roughly accounted for by the last layer of optical thickness of 0(/3) 
(equivalently G(£p) thickness in the physical metric). 

3.4 The optical expander 

We will now make a basic assumption regarding the structure of the A-stretched 
horizon. We declare that the geometrical structure of X is still relevant to characterize 
the interactions of strongly-coupled degrees of freedom building up the stretched horizon. 
More precisely, we model the stretched horizon as an effective hyperbolic lattice of O{0) 
spacing with a worth of 0(1) Q-bit degrees of freedom per lattice point. In other words, 
we have a strongly-coupled discrete quantum system on a cut-off expander graph. 



The discussion leading to ( 20 ) generalizes to any non-trivial CFT perturbed by relevant and irrelevant 



operators with respective scales M and A. 
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The dimension of this Hilbert space is fixed by the matching to the Bekenstein- 



Hawking entropy. In particular, we have the extensive behavior (24), now applying to the 



dimension of the Hilbert space, rather than simply the entropy of excitations: 

logdimHsH ~ V^sh ~ 5 BH , (27) 



where Vsh ~ K P and we have used (26). 



Setting V = V~ = (3 d , we deal with a single thermal cell, and the associated entropy 



is just the specific entropy S e s = S ce xi = N. Inverting (25) we find the relation between 



the maximal depth of the stretched horizon and the specific entropy 

= log N , (28) 

for the total optical depth at the 'bottom' of the stretched horizon. 

At this point, we can make contact with the discussion in section 2, characterizing 
horizon scrambling in terms of a discrete diffusion process with step /3, playing on the 
expander graph cutoff to depths z\ < z < zp. The expander-diffusion process scrambles 
0(n) degrees of freedom of a single thermal cell in a time of order /31og(n), as the random 
walk reaches the 'bottom' of the stretched horizon at z = zp. Any subsequent scrambling 
over scales larger than a thermal cell of the 7^ metric proceeds by standard diffusion on 
the 'bottom' of the expander with induced (optical) metric 



ft 

zp 



^ ^ T y dyW = N 2/d Tij dyW . (29) 



Hence, a length scale L 3> (3 in the 7^ metric corresponds to a length scale L N 1 ^ 
{NV) l l d in the bottom metric (29). Such a patch is slow-scrambled in a time of order 



f3~ l (VN) 2 l d . The number of thermal cells on the patch is V/f3 d . Since iV is the entropy 
per thermal cell we may write the total entropy of the patch as S ~ iV V f3~ d , so that the 
total slow-scrambling time is the standard diffusive /3 S 2 ^ d result, as in It is important 
to notice that this result is non-trivial, since it has been derived for systems which have 
O(N) degrees of freedom per thermal cell, rather than 0(1) degrees of freedom as was 
the case in the Introduction. 

Combining these estimates we finally find the fast-scrambling law 

t s (L) ~ /3 log N , L < (3 (30) 
below the thermal cell scale, and the slow-scrambling law 

r s (L)~(3S 2/d , L>/3 (31) 
for distances well above the thermal cell. 



The estimate (31 ) is quite general, since it only depends on the bottom density of the 



expander. On the other hand, (30) assumes tacitly that the expander graph starts at 
z = za ~ 0, i.e. the stretched horizon covers the whole Rindler region. Such a situation 
is only appropriate for highly stringy black holes dual to weakly-coupled CFTs. Smooth 
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one 7— thermal cell 



Figure 2: Schematic view of the effective kinetic scrambling model in the optical frame, fea- 
turing a Cayley tree picture of the stretched horizon, a hyperbolic 'Rindler atmosphere' and 
an asymptotic AdS region of optical depth of 0(/3). The total optical depth of the hyperbolic 
section [0, zp] is O (f3 log N). A localized classical probe injected from the asymptotic boundary 
( vertical blue line) is either absorbed or reflected at z = z\. If reflected (curved red line), it 
scrambles fast in the effective chaotic billiard of the Rindler atmosphere. If absorbed (jagged 
red line), it scrambles fast by diffusion on the expander graph until it reaches the bottom at zp. 
Once at the bottom, it scrambles slowly by diffusion on the bottom. The fast-scrambling patch 
measured by the 7^ metric is always about one thermal cell. 
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horizons with a large Rindler region where LEFT is valid support 'thin' stretched horizons 
for which the fast-scrambling diffusion operates for a smaller time of order /31og(iV) — z\. 
In particular, in the limit in which we push the LEFT validity down to Planckian distances 
from the horizon, we have za ~ z p up to 0(1) factors, and the optical thickness of the 
stretched horizon is only of 0((3). In such a situation the random walk diffusion is always 
slow-scrambling. However, as mentioned above and argued in [7] , surface-scattering effects 
turn the LEFT region < z < zp into a chaotic billiard with Lyapunov time of order 
(3 log N. Therefore, the combined effect of classical chaos with 0(1) 'atmospheric' collision 
events and the standard diffusion upon absorption by the expander graph account for the 
expected scrambling properties of idealized probes. A cartoon of this picture is offered in 
figure 2. 

4 Beyond the probe approximation 

Up to this point have modeled scrambling in a rather naive fashion, using the random 
walk model. Our main contention is the ability to perform ballistic-fast scrambling in this 
way, provided we define an appropriate background geometry for the diffusion process. 
This picture requires the operational definition of a probe system which is stable along 
the time evolution of the system and effectively classical on a coarse-graining time scale 
(3. It would be interesting to free the discussion from this requirement, in such a way that 
the fast-scrambling property is characterized in a more fundamental way. 

In this section we advance some preliminary remarks about thermalization in quantum 
systems defined on expander graphs. We make no explicit reference to a probe subsystem 
or any other quasiparticle degree of freedom. Rather, we focus on generic properties of 
entanglement with respect to bi-partitions of the expander graph, and point out some 
difficulties in arriving at a good physical definition of scrambling. 

Let G be a connected graph where a local quantum system is defined and A U A a 
bi-partition into two connected components, with the convention that A has the smaller 
vertex size: N V (A) < N/2. Given a general state on G, described by a density matrix p, 
we define the marginal density matrix of A and its entanglement entropy as 

p A = Tta P , Sa = -Tr^ Pa log pa ■ 

A natural criterion to declare a system as scrambled is to demand that Sa be bounded 
below by some fixed extensive function, for all choices of A under the stated conditions, 

S A > C N V (A) , N V (A) < N/2 . (32) 

For systems with finite-dimensional Hilbert space at each vertex, this criterion is equiva- 
lent to maximal mixing of pa when the constant C saturates to its upper bound C max = 
log dim r H, v . In other cases, it incorporates the fact that we may have different effective 
temperatures depending on the total energy of the state. Alternatively, one can speak of 
maximal mixing on the fixed energy surface, assuming ergodicity, i.e. no other conserved 
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quantities besides the total energy. The intuition behind the criterion is that proper sub- 



systems see the complement as a thermal reservoir. The condition (32) is expected to 
be generic [22] for randomly chosen global states p, and we shall refer to it as the Page 
criterion. 

The scrambling time may be defined as the typical lapse of Hamiltonian evolution 
needed to satisfy the Page criterion, starting from an initial state which explicitly violates 
it. Of course, most of the difficulties in establishing general results lie on the lack of 
precision in just every concept appearing in the previous sentence. 

In order to sharpen the language a bit we may focus on globally pure states on G. 
Then, a good choice for the 'least thermalized' initial state may be a completely factorized 
one ® v \ip v ) V) for which Sa = for all partitions. An up-to-date discussion of scrambling 
times for such initial states can be found in [T3] . 

Another natural choice of initial state is one whose entanglement properties imitate 
those of the ground state. For systems with short-range correlations, such ground states 
satisfy the so-called 'area law' for the entanglement entropy, i.e. Sa ~ size[cM], where 
size[cM] can be measured in terms of the edge size of the boundary, Nt( OA), for fc-regular 
graphs. The usual intuition about local systems in Euclidean geometry suggests then that 
the scrambling consists on the gradual evolution from 'area law' to 'volume law' for the 
entanglement entropy. 

Already at the level of these primitive notions, we see that expander graphs are special 
in their behavior with respect to extensivity. In particular, there is no clear-cut distinction 
between 'area-law' and 'volume-law' behaviors when G is an expander. As a consequence 



of the expansion property (14), any area-law entanglement entropy of a bi-partition G 



A U A satisfies the Page criterion 

S A = s A N t ( OA) > h k s A N v ( A) 

from the beginning. Such 'area- law' states are not difficult to engineer and do not look 
'thermal' at all. To illustrate this fact, let us consider a simple Hamiltonian with localized 
vertex terms and 'hopping' interactions for each edge of the graph, 

H= h *+ £ h i> ( 33 ) 

u£ vertices edges 

where h v acts on the single-vertex Hilbert space H v and the hopping term hg acts on the 
two- vertex Hilbert space %t = % v ® % w associated to the edge £ = (vw). Let us suppose 
now that hi <C h v and construct the ground state in the hopping expansion in powers 
of A, a small parameter characterizing the matrix elements of hi. To zeroth order, the 
vacuum is determined by the on-vertex vacua |0)„ for each operator h v , 

| Vac) = | n ) + 0(A) = ® v \0) v + 0(A) , 

a completely factorized state for which Sa = for any bi-partition. To leading order in 
the hopping terms the ground state has the form 

|Vac) = a|n) + A ^ K) + 0(A 2 ) 

l £ edges 
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where the sum extends over all links i = (vw) of the lattice G and the state \£) differs from 
the completely factorized on-site vacuum |IT) by an excitation on the £-edge Hilbert space 
Hi = l-i v ®l-L w . It is orthogonal to the zeroth-order vacuum, ( II 1 1 ) = 0, and its particular 
form is determined by the matrix elements of hi through the usual formulae of first-order 
perturbation theory. Normalization sets the constant a to the value a = \J\ — X 2 Ni, with 
Ng the total number of edges in G. 

Given now a bi-partition G = A U A, we define the state 

\\ A o) = ^\n A o) + Ay^K), 

where the link sum includes only those edges lying entirely within the interior A° of the 
subgraph A and | H A °) = ®veA°\tyv Using this state we can write the vacuum to leading 
order in the hopping expansion as 

| Vac) = \X A o)(g)\X A o) + X^2\e a ) + 0(X 2 ) , 

ta 

where £q denotes edges on the boundary OA. In other words, the vacuum remains factor- 
ized except for 0(A) non-trivial excitations localized on the boundary of the partition. It 
follows that the entanglement entropy of the partition is proportional to the number of 
edges on the boundary, i.e. 

S A = XcN e {dA) + 0{X 2 ) , 

with c a positive numerical constant. 

Hence, we recover the expected result: ground states with very short-range correlations 
have 'area-law' entanglement across bi-partitions. The peculiar behavior of expander 
graphs regarding Page-like criteria can be further emphasized by looking at the opposite 
limit, in which we regard the vertex hamiltonian density as negligible compared to the 
hopping term, i.e. in the limit \h v \ <C \hi\. Let us assume that the Hamiltonian is 
homogeneous, in the sense that hi is given by the same operator for all links, so that 
[hi, hgi] = for all links on G. Then, the evolution operator can be approximated as 

U G (t)= J] U *> U e = e- m *, (34) 

l £ edges 

with hg independent of £, so that the ordering in the operator product is actually immate- 
rial. If the initial state is still given by the factorized ground state of the h v Hamiltonian, 
S8)„|0) v , since in general [h v , hi] ^ when v G £, we see that each link operator U vw will 
generate entanglement in the % v ® % w subspace in a time scale set by the eigenvalues of 
hi. Since hi is interpreted as the hopping operator, we call j3 the characteristic time scale 
for its inverse normal frequencies. The resulting state with 0(1) entanglement across 
links is generally termed a graph state in the literature (cf. |23j), and its entanglement 
entropy is of the order of the rank of the off-diagonal adjacency sub-matrix connecting A 
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and A. For a regular graph with fixed coordination, k, this so-called 'Schmidt rank' is of 
the order of the boundary edge-size, Nt(dA), up to fc-dependent degeneracy factors. 

We thus conclude that hopping-dominated hamiltonians on regular graphs achieve the 
graph state with area-law entanglement in times of 0(/3). If the regular graph is also an 
expander, then we satisfy the Page criterion in 'scrambling' times of 0(/3). Such a ultra- 
fast scrambling time is undoubtedly a consequence of the highly integrable properties of 
the pure-hopping Hamiltonian with commuting link operators. For example, taking % v 
to be a two-level quantum system and hi determined by a single fundamental frequency 
of time scale 0, the evolution operator (34) is actually periodic on time scales of (9(2/3), 
so that it has very small times of Poincare return, making it hard to talk physically about 
'scrambling' or any type of thermalization for that matter. Q 

It would be interesting to study the scrambling times, in the sense of the Page cri- 
terion, for more generic Hamiltonians in the class (33), without the extreme simplifying 
assumptions of vertex and/or link domination. In any case, the comments offered in this 
section suggests that the Page criterion is not nearly enough to successfully characterize 
scrambling times in expander graphs. 



4.1 Going ballistic 

Given the previous comments, it might be interesting to look for alternative charac- 
terizations of 'fast' scrambling in expander graphs. One possibility is to focus attention 
on the 'ballistic' character of expander-diffusion as the defining property, and look for 
quantum criteria adapted to this fact. The diffusion model based on random walks as- 
sumes the stability of the walking probe. This in turn is equivalent to the assumption that 
the system has a conserved quasiparticle density or, more generally, a locally conserved 
current. On general grounds, retarded correlation functions of such local currents will 
show a characteristic 'diffusion pole' in the lower half plane of complexified frequencies. 
For local diffusion on an Euclidean geometry with diffusion constant D, such a pole takes 
the form u^is = —i D k 2 + . . . in a large wavelength expansion. Setting \ k\ ~ 1/L for 
a fluctuation of wavelength L, we find the standard diffusion time scale D~ l L 2 of a slow 
scrambler J3 

According to the general discussion of section 2, the behavior of retarded correlation 
functions is expected to be quite different when the diffusion takes place in an expander 
graph. As a preliminary discussion, we may look at the continuum approximation and 
concentrate on the behavior of retarded correlation functions for field theories on hy- 
perbolic manifolds. A first indication in this direction is the behavior of the free-field 
conformal Green's functions on hyperbolic space-times of the form R x H n , with the first 

7 The same can be said about the random Ising model example presented in [13]. In that case the 
scrambling time is of 0(log N) because it takes 0(log N) time steps to build a random lattice with 
0(log N) average vertex degree. 

8 This diffusion quasi normal mode is visible in the gravitational description of hydrodynamics, even 
if a quasiparticle picture is absent in this framework, provided D is of 0(1) in the large N limit. 
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factor representing a time-like dimension]^] 



n — 1 



G ret (t } x) = 0, # - + — — i? H « t, x ) ~ 5(t - L x 



1 



,n-2 ' 



4n 7 / u " (sinhL-, 

/ret 

(35) 

where (t, x) G R x H n and is the hyperbolic geodesic distance from the conventionally 
chosen origin of coordinates on H n to the point x (cf. for example [IS]). We see that 
the retarded correlation function is indeed ballistic, since it is strictly supported over the 
past light-cone of the origin, despite the presence of a constant-curvature term which 
simulates a mass. This Green function is automatically thermal, since it is periodic under 
imaginary-time shifts t — > t+ 2m when integrated against arbitrary sources. The effective 
temperature for a hyperboloid of radius R is then (3~ l = (2irR)~ 1 . 

The fact that perturbations on a free system behave ballistically at finite tempera- 
ture is not surprising. After all, the retarded Green function of any free field theory is 
a c-number, being proportional to the commutator function of free fields, and thus it is 



the same in any state. On the other hand, the functional form in (35) is largely dic- 
tated by conformal symmetry, and in particular the hyperbolic Green function must be 
proportional to the vacuum Green function on Minkowski space, since R x H™ can be 
conformally mapped into a proper region of R x R ra , given by the causal development of 
the unit ball B n on the t = spatial section R n . This region of Minkowski space, denoted 
X>(B n ), has the form of a 'diamond' based on B n . 



This fact suggests a very interesting and non-trivial generalization of (35) to the case 
of fully interacting (but still conformal) field theories. It is known (cf. for example [24] ) 
that a CFT thermal state of temperature 1/2tv on R x H ra is conformally related to a 
vacuum state on Rx R™ in the following sensej^] any time-dependent correlation function 
on the hot hyperboloid is conformally related to a vacuum Minkowski correlation function, 
restricted to T>{Q n ). 

Denoting rj^ the Minkowski metric on R x R n and the hyperbolic metric on 
R x H™, we write rj^ = Q 2 for the conformal rescaling. It turns out that the whole 
hyperbolic space-time is mapped to the finite domain T>{Q n ) of Minkowski space-time. 
The Minkowski vacuum, as measured by local correlation functions restricted to D(B n ), 
is equivalent to the mixed state obtained by tracing over the degrees of freedom 'outside' 
the unit ball, i.e. we may replace the Minkowski vacuum by the entanglement density 
matrix on the unit ball. The conformal map sends this entanglement density matrix into 
the thermal density matrix of the hyperbolic space-time. If the correlation functions are 
defined for scaling operators of conformal dimension Aj we have 

( dfoO ■ ■ ■ O s {y s ) ) j8>H n = ft(yi) Al • • • n(y s ) A ° ( O x (y x ) ■ ■ ■ O s (y s ) ) vac , c(B ") (36) 

where the operators 0(y) at points y G R x H n are the images under the conformal map 
of the operators 0(y) at points y G V(B n ). 



9 When comparing formulae with previous sections, notice that we have set n = d+1 in this subsection. 
10 Notice that a variant of the same argument was used before to relate the Hartle-Hawking state on 
Rindler to the thermal state on the hyperboloid. 
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Consider now the retarded two-point function of an operator 0(y) = 6&(f,x) in the 
Minkowski vacuum of the CFT. We may write it as the imaginary part of the time-ordered 
correlation function, i.e. 

G vct (i,x) oc 9(-t) • Im <T [<5 A (0,0) & A (i, x)]) vacV = Im ( __ 2 } + - Q)A , 

where 9(t) is the Heaviside function and we have used conformal invariance in the last 
step. Since we assume the CFT to be unitary, A is real, and the retarded correlation 
function is a distribution supported on the Minkowski (past) light cone of the origin, 



t 2 = \x\ 2 , t < 0, just like it was for free fields. Applying now the conformal map (36) we 
learn that the thermal retarded Green function on the hyperbolic manifold 

G Iet {t,x) = n(t,x) A fi(0,0) A G Iet {i,x) 

is supported on the image of the light-cone under the conformal map. The conformal 
transformations send light-cones into light-cones, so that G rct (t,x) is supported on the 
past light-cone of the origin in R x H n . This proves that the fully interacting Green 
function in the hyperbolic space-time is supported on the light-cone, and thus any diffusion 
poles (for A (y) equal to a conserved current of dimension A = n) must be compatible 
with the ballistic propagation on the hyperboloid. 

It would be interesting to prove the ballistic character of thermal Green functions 
directly for discrete quantum systems on finite expander graphs. Perhaps such a proof 
can be approached by first solving the intermediate case of the Bethe lattice, a discrete 
expanding geometry of infinite volume. 



5 Concluding remarks 

In this paper we have pointed out that a classic result in graph theory, namely the 
exponentially fast mixing of random walks on expander graphs, provides an interesting 
semi-phenomenological model for fast scrambling at black-hole horizons. 

Let us briefly recapitulate the logical steps leading to this conclusion. First, expander 
graphs are discrete versions of hyperbolic geometry, so that fast brownian mixing is akin 
to fast diffusion on hyperbolic spaces. On the other hand, regular horizons are confor- 
mal to asymptotically locally hyperbolic spaces, the Hartle-Hawking quantum state being 
mapped to a homogeneous thermal ensemble on this so-called 'optical frame'. Working 
in this frame, it becomes then natural to identify the stretched horizon with a cut-off 
expander graph, hosting a quantum system with effective temperature 1//3. This identifi- 
cation relates fast-scrambling to diffusion in the 'optical depth' of the graph, turning into 
more standard 'slow' scrambling when the bottom of the expander is reached. The rele- 
vant length scale separating the two scrambling modes is naturally given by the thermal 
length of the system. 

This geometrical picture of the stretched horizon also dovetails well with the chaotic- 
billiard model of 'atmospheric' scattering developed in [7]. While the expander graph 
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diffusion accounts for the scrambling of information absorbed by the stretched horizon, the 
billiard model accounts for the scrambling of information kept in the thermal atmosphere. 

The most important open problem regarding these ideas is the derivation of the ex- 
pander graph structure from a complete holographic description. In this regard, a natural 
expectation might be that the expander graph structure could survive the limit to weak 
coupling in the holographic dual, since the fast-scrambling scale, proportional to log N, is 
not explicitly dependent on 't Hooft coupling factors. If this expectation is borne out, the 
fast-scrambling properties of sufficiently non-trivial matrix models (cf. [23]) could rely 
on identifying an expander graph structure in its Hilbert space, perhaps allowing for the 
expander graph to be defined SIS db random lattice. 

For this purpose, it is important to free the formal treatment of scrambling from the 
simple random-walk model. A preliminary discussion along those lines is offered in the 
fourth section of this paper, where it is argued that standard formal criteria, such as the 
Page test, are not necessarily sharp enough, and a more formal characterization in terms 
of ballistic behavior of thermal correlation functions is proposed. 

It should be emphasized that the local operators referred to in section 4 are bulk 
objects, with an extremely non-local relationship to a given basis of local operators on the 
holographic dual. Indeed, the formulation of scrambling in the 'matrix-model', language 
brings into focus the question of what observables are adequate. On general grounds, 
correlation functions of local operators in large- iV CFTs are expected to show quasi normal 
behavior with decay time of 0(/3) and long-time quantum noise with infinite-time average 
of 0(e~ s ) and associated Poincare recurrences on time scales of 0(/3e s ) (cf. [26]). 
Quasinormal behavior is thus irrelevant for time scales longer than 0(/3S), becoming then 
subsumed into the quantum noise. On the other hand, diffusion times of slow-scramblers 
can be detected in the correlation functions of locally conserved densities, associated to 
the standard diffusion poles in any hydrodynamical description. 

Detecting the log N time scale directly in CFT correlation functions is likely to require 
a non-trivial CFT definition of local operators, deep inside the near-horizon region. How 
to characterize such extremely bulk-local operators remains a difficult challenge. Other 
evidences for fast thermalization times are directly visible at iV = oo in terms of non-local 
sharp probes, as in [27]. 
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